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Abstract.  The  occurrence  rate  of  great  geomagnetic 
storms  displays  a  pronounced  semiannual  variation.  Of 
the  forty-two  great  storms  during  the  period  1940-1990, 
none  occurred  during  the  solstitial  months  of  June  and 
December,  and  40%  (17)  occurred  during  the  equinoctial 
months  of  March  and  September.  This  suggests  that  the 
semiannual  variation  found  by  averaging  indices  is  not  the 
result  of  some  statistical  effect  superposed  on  the  effects 
of  random  storm  occurrence  but  rather  is  dominated 
by  the  storms  themselves.  Recent  results  indicate  that 
the  intense  southward  interplanetary  magnetic  fields 
(IMFs)  responsible  for  great  storms  can  reside  in  the 
postshock  plasma  preceding  the  driver  gas  of  coronal 
mass  ejections  (CMEs)  as  well  as  in  the  driver  gas 
itself.  Here  we  propose  that  strong  southward  fields  in 
the  postshock  flow  result  from  a  major  increase  in  the 
Russell-McPherron  polarity  effect  through  a  systematic 
pattern  of  compression  and  draping  within  the  ecliptic 
plane.  Differential  compression  at  the  shock  increases 
the  Parker  spiral  angle  and,  consequently,  the  ausimuthal 
field  component  that  projects  as  a  southward  component 
onto  Earth’s  dipole  axis.  The  resulting  prediction  is 
that  southward  fields  in  the  postshock  plasma  maximize 
at  the  spring  (fall)  equinox  in  CMEs  emerging  from 
toward  (away)  sectors.  This  pattern  produces  a  strong 
semiannual  variation  in  postshock  IMF  orientation  and 
may  account  at  least  in  part  for  the  observed  semiannual 
variation  of  the  occurrence  of  great  geomagnetic  storms. 

Introduction 

Geomagnetic  storms,  defined  as  periods  of  strengthened 
ring  current  in  the  inner  magnetosphere,  correlate 
remarkably  well  with  southward  excursions  of  the  inter¬ 
planetary  magnetic  field  (IMF)  [e.g.,  Burton  et  al.,  1975]. 
Thus  predicting  storms  becomes  a  matter  of  predicting 
IMF  orientation.  In  general  storms  arc  attributed  to  two 
heliospheric  features  recurrent  stream  interaction  regions 
and  transient  coronal  mas.=  ejections  (C.MEs).  The  latter, 
when  ejected  fast  enough  to  form  a  shock  wave,  produce 
the  largest  storms.  For  example,  Gosling  et  al.  [1991]  find 
that  13  of  the  14  largest  storms  as  measured  by  Kj,  were 
caused  by  CMEs  with  shocks.  Tsurutani  et  al.  [1988, 
1992]  demonstrate  that  the  southward  fields  that  cire  the 
direct  cause  of  CME-driven  storms  reside  cither  in  the 
postshock  plasma  that  precedes  the  driver  gas  or  in  the 
driver  gas  itself  or  both.  To  predict  the  largest  storms, 
then,  one  wishes  to  predict  the  orientation  of  the  IMF  in 
the  postshock  plasma  and  driver  gas  of  fa.st  CMEs. 

A  clue  to  predicting  IMF  orientation  in  po.stshock 
plasma  and  driver  gas  comes  from  the  semiannual 
variation  of  large  storncs.  The  next  section  documents 
the  remarkable  strength  of  this  variation.  The  following 
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section  demonstrates  how  the  semiannual  variation  can 
arise  from  the  heliospheric  ordering  of  postshock  flow  and 
gives  a  scheme  for  predicting  IMF  orientation  there.  The 
final  section  discusses  other  possible  contributions  to  the 
semiannual  variation,  including  IMF  orientation  in  the 
driver  gas. 

Semiannual  Variation  of  Large  Geomagnetic  Storms 


Xo 


The  semiannual  variation  of  geomagnetic  activity  is 
well  known  and  has  been  analyzed  by  a  number  of  meth¬ 
ods  (see,  e.g.,  Crooker  and  Siscoe  [1986a]  or  Gonzalez  et 
al.  [1992]  for  review) .  It  is  usually  treated  as  a  statistical 
effect  and  attributed  to  a  mechanism  that  gives  stronger 
solar  wind  -  magnetosphere  coupling,  on  the  average,  in 
spring  and  fall.  For  example,  Russell  and  MePherron 
[1973]  attribute  the  semiannual  variation  to  the  fact 
that  the  average  IMF,  oriented  at  the  Parker  spiral 
angle  in  the  solar  equatorial  plane,  projects  a  southward 
component  at  Earth  in  the  coupling-effective  geocentric 
solar  magnetospheric  coordinate  system  whenever  its 
polarity  points  toward  the  sun  in  spring  and  away  from 
the  sun  in  fall.  The  southward  component  achieved  by 
this  projection  at  equinox  ranges  from  28%  to  40%  of  the 
IMF  field  strength,  depending  upon  universal  time. 

Less  known  is  the  fact  that  the  occurrence  rate  of  large 
geomagnetic  storms  also  shows  a  semiannual  variation, 
which  cannot  be  explained  by  the  effects  of  average  fields. 
Newton  [1948]  points  out  that  the  monthly  distribution 
of  great  storms  at  Greenwich-Abinger  during  the  period 
1875-1946  has  marked  peaks  in  spring  and  fall.  Figure  la 
shows  a  histogram  of  the  data  used  by  Newton,  updated 
to  1954  [Jones,  1955].  Of  the  total  112  great  storms,  30% 
occurred  in  the  peak  equinoctial  months  of  March  and 
August  and  6%  during  the  solstitial  months  of  June  and 
December.  Even  the  sub-set  of  the  11  most  outstanding 
storms  shows  a  semiannual  variation.  Allen  [1982]  obtains 
a  similar  result  by  plotting  the  monthly  distribution  of 
major  .storms  determined  by  the  parameter  A’  >  80 
during  the  period  1932-1980  [sec,  also,  Gonzalez  et  al., 
1992].  Figure  lb  updates  his  result  and  restricts  it  to 
larger  storms,  with  A*  >  100,  during  the  period  1932- 
1989.  Of  the  total  140  storms,  33%  occinrcd  in  the  peak 
equinoctial  months  of  March  and  September,  and  3.5%  in 
the  solstitial  months  of  June  and  December  (all  in  June). 

Finally,  Figiuc  Ic  gives  the  monthly  distribution  of 
great  storms  defined  by  the  limit  £>,,(  <  -250  nT, 
following  Tsurutani  et  al.  [1992].  Since  D.,i  dates  back 
only  to  1957,  we  extend  the  list  of  storms  to  include  the 
period  from  1932  by  using  the  index  A*  as  a  proxy  for 
Dsi  prior  to  1957  in  the  following  manner.  For  storms 
with  D31  <  —250  iiT  during  1957-1990,  72%  (18/25)  have 
A*  >  160.  Conver.sely,  of  storms  with  A*  >  160  during 
this  period,  69%  (18/26)  have  D,i  <  -250  nT.  Thus  we 
choose  A*  >  160  as  a  reasonable  measure  of  great  D,i 
storms  for  years  prior  to  1957  and  lengthen  the  list  from 
25  to  42  cases.  Of  these,  40%  occurred  during  the  peak 
equinoctial  months  of  March  and  September,  and  none 
during  the  solstitial  months  of  June  and  December. 

The  pronounced  semiannual  variation  in  the  histograms 
in  Figure  1  docs  not  arise  as  the  result  of  large  storms 
occurring  randomly  from  one  year  to  the  next  in  one 
or  the  other  of  the  equinoctial  periods.  Rather,  the 
variation  is  remarkably  well-defined  even  within  single 
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Fig.  1.  The  semiannual  variation  of  the  occurrence  of 
great  geomagnetic  storms  a.)  with  D  range  >  60'  or  H  or 
Z  range  >  300  nT  at  Greenwich  or  Abinger  [Jones,  1955), 
b.)  as  measured  by  A*  >  100  nT,  and  c.)  by  D^t  <  -250 
nT,  after  Tsurutani  et  al.  [1992].  The  indicates  that 
A’  >  160  nT  was  used  as  a  proxy  for  Du  prior  to  1957. 

years  that  have  large  storms.  Figure  2  demonstrates 
this  fact  with  plots  of  the  variation  of  monthly  averages 
of  D,t  for  those  years  when  any  single  monthly  average 
is  <  -50  nT.  In  all  three  cases,  there  arc  pronounced 
peaks  near  both  equinoxes.  The  phase  varies  from  year 


Fig.  2.  The  semiannual  variation  of  monthly  averages  of 
D,t,  with  scale  inverted,  for  those  years  beginning  with 
1957  in  which  any  monthly  average  was  more  negative 
than  -50  nT. 


to  year-for  example,  the  .spring  peak  occurs  in  February, 
March,  and  April,  in  the  years  1982,  1957,  and  1960, 
respectively -but  the  overall  clarity  of  the  variation  is 
preserved.  Bartels  [1963]  finds  similar  patterns  for  the 
indices  Kp  and  Aj,  during  disturbed  years.  The  variable 
phase,  also  founcl  in  the  power  spectral  analysis  of  Ap 
(Gonzalez  et  al.,  1992],  may  reflect  the  remnant  of  CME 
occurrence  variability  after  modulation  by  the  mechanism 
that  creates  the  semiannual  wav .. 

The  results  in  Figures  1  and  2  demonstrate  semiannual 
variations  with  amplitudes  that  far  exceed  those  obtained 
by  averaging  magnetic  indices  over  many  years.  For 
example,  Legrand  and  Simon  [1989]  point  out  that 
the  semiannual  variation  in  the  aa  index  averaged  over 
the  period  1868-1980  hcis  an  amplitude  equal  to  13.5% 
of  the  average  activity  level,  whereas  the  occurrence 
ft’equency  of  severe  storms,  defined  by  Aa  (daily  sum 
of  aa)  >  100  nT,  shows  a  modulation  of  80-90%.  This 
amplitude  difference,  as  well  as  the  fact  that  individual 
yeairs  that  have  great  storms  show  the  most  pronounced 
semiannual  variation,  as  in  Figure  2,  imply  that  storms 
are  responsible  for  most  of  the  semiannual  variation, 
rather  than  its  being  a  statistical  effect  that  one  can 
sift  out  of  the  data  by  averaging  over  the  effects  of 
random  storm  occurrence.  Green  [1984]  reaches  the  same 
conclusion  based  on  his  analysis  of  monthly  averages  of 
aa,  designated  Aq.  He  shows  that  the  amplitude  of 
the  semiannual  variation  increases  with  increasing  Aq 
and  that  its  modulation  in  time  is  nearly  identical  to 
the  time  modulation  of  the  amplitude  of  the  semiannual 
variation  in  Aq  occurrence  for  stormy  days  (Aq  >  30). 
Although  the  semiannual  variation  is  apparent  at  all  levels 
of  activity  [e.g.,  Bartels,  1963;  Green,  1984;  Gonzalez  et 
al.,  1992],  its  dominance  at  high  levels  clearly  indicates 
that  some  mechanism  operates  to  occasionally  produce 
very  strong  southward  IMF  at  equinox  and  not  at  solstice. 

Postshock  Magnetic  Field 

Figure  3  presents  a  schematic  drawing  of  the  cross- 
section  of  a  coronal  mass  ejection  (CME)  deforming  the 
spiral  magnetic  field  pattern  in  the  ecliptic  plane  as  it 
moves  outward  from  the  sun  faster  than  the  ambient 
medium.  The  CME  presumably  emerged  from  a  helmet 
streamer  [e.g.,  Hundhausen,  1988]  and  expands  into  the 
solar  wind  along  a  sector  boundary,  marked  by  the 
heliospheric  current  sheet  (HCS)  dividing  fields  pointing 
away  from  and  toward  the  sun.  Because  of  the  spiral  field 


Fig.  3.  Cross-section  in  the  ecliptic  plane  of  a  coronal 
mass  ejection  (CME)  expanding  into  the  heliosphere 
along  a  sector  boundary,  marked  by  the  heliospheric 
current  sheet  (HCS).  Spiral  magnetic  field  lines  from  the 
leading  sector  bend  at  the  shock  and  drape  around  the 
CME,  effecting  an  increase  in  the  Parker  spiral  angle  and 
the  field  stren^h  in  the  postshock  region. 
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geometry,  the  leading  edge  of  the  CME  becomes  mostly 
that  portion  west  of  its  intersection  with  the  HCS  (see, 
also,  Gosling  et  al.  [1987]).  The  shock  wave,  represented 
by  a  dotted  line,  compresses  the  component  of  the  spiral 
field  that  lies  perpendicular  to  the  shock  normal.  Thus 
it  increases  the  Parker  spiral  angle  as  well  £is  the  total 
strength  of  the  field  from  the  leading  sector  across  most  of 
the  leading  edge  of  the  CME.  Further  field  intensification 
takes  place  from  the  shock  to  the  CME  boundary  as  the 
field  stretches  around  the  oncoming  CME,  in  analogy  with 
magnetic  field  behavior  in  Earth’s  magnetosheath  [e.g., 
Crooker  et  al.,  1982]. 

Figure  4  shows  the  IMF  draped  over  the  CME  viewed 
toward  the  sun.  The  pattern  is  taken  from  the  model 
calculated  for  Earth’s  magnetosheath  with  an  upstream 
IMF  directed  30°  from  the  Earth-sun  line  [Crooker  et 
al..  1985].  (A  pattern  for  the  45°  Parker  spiral  angle 
at  Earth  is  not  available.)  The  length  of  each  vector  is 
proportional  to  field  strength.  The  arrowheads  have  been 
omitted  so  that  one  can  use  the  diagram  to  represent 
either  toward  or  away  polarity.  Earth’s  position  could 
be  anywhere  relative  to  the  CME  but  is  shown  as  a 
dot  on  the  CME’s  ecliptic  equator  near  the  region  of 
maximum  field  strength.  The  draped  field  there  lies  in 
the  geocentric  solar  ecliptic  (GSE)  X-Y  plane,  primarily 
in  the  Ycse  direction  (parallel  to  the  solar  equatorial 
Y  direction  at  equinox').  Earth’s  Z  axis  in  geocentric 
solar  rnagnetospheric  (GSM)  coordinates  is  tilted  by 
the  maximum  value  of  35°  in  the  sense  appropriate  for 
September  equinox  at  1030  UT.  For  an  IMF  directed 
away  from  the  sun,  as  in  Figure  3,  and  as  indicated  by 
the  large  horizontal  vector  in  Figure  4,  its  projection 
onto  the  Zgsm  axis  provides  a  southward  component 
equal  to  57%  of  the  Ycse  component,  which,  according 
to  Figure  3,  is  nearly  the  entire  IMF  vector.  Thus  the 
Russell-McPhcrron  [1973]  effect,  applied  to  the  geometry 
of  the  postshock  field  preceding  a  CME,  provides  a  much 
stronger  southward  IMF  component  than  in  its  original 
application  to  the  average  Parker  spiral  field.  FVom  it  one 
can  predict  that  postshock  plasma  will  be  gcoeffcctive 
in  northern  hemisphere  spring  (fall)  upon  approach 
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Fig.  4.  View  toward  the  sun  of  an  array  of  magnetic 
field  vectors  draped  over  an  approaching  CME,  from  the 
magnetasheath  draping  model  of  Crooker  et  al.  [1985]. 
Superposed  on  the  array  arc  vectors  illustrating  the 
{xjstshock  Russcll-McPhcrron  effect.  From  Earth,  marked 
by  a  dot,  extends  a  vector  pointing  in  the  Ycse  direction 
representing  the  compressed,  draped  field  preceding  the 
CME  in  Figiire  3.  Its  projection  onto  the  ZesM  a^xis, 
shown  tilted  the  maximum  of  35°  from  the  GSE  system 
in  the  sense  appropriate  for  September  equinox,  forms  a 
geocffcctive  southward  component. 


to  toward-to-away  (away-to-toward)  sector  boundaries, 
assuming  CMEs  emerge  from  helmet  streamers  there,  as 
in  Figure  3.  Without  this  assumption,  the  prediction  still 
holds  for  CMEs  from  toward  sectors  in  spring  and  away 
sectors  in  fall. 

It  is  important  to  note  that  the  draping  invoked  here  to 
produce  southward  IMF  concerns  only  those  fields  lying  in 
the  ecliptic  plane,  in  contrast  to  the  proposal  by  Gosling 
and  McComas  [1987]  that  draping  lifts  the  ambient  field 
out  of  the  ecliptic  to  produce  southward  IMF.  Models 
of  Earth’s  magnetosheath  field  indicate  that  substantial 
out-of-the-ecliptic  components  arise  only  when  the  IMF 
is  nearly  radial  [Kartelev  and  Mastikov,  1982;  Crooker  et 
al.,  1985].  For  exactly  radial  IMF,  the  field  diverges  from 
the  nose  of  the  magnetosphere,  or  CME,  in  this  Ccise,  and 
drapes  directly  upward  and  downward  in  the  meridian 
plane  containing  the  nose,  £is  pictured  by  Gosling  and 
McComas.  However,  for  the  IMF  spiral  angle  of  30°  used 
to  construct  Figure  4,  the  divergence  point  in  the  field 
pattern  against  the  CME  surface  has  moved  off  to  the  left, 
leaving  draped  field  vectors  lying  mostly  parallel  to  the 
ecliptic  plane,  especially  in  the  region  of  maximum  field 
strength.  Since  at  Earth’s  orbit  the  ambient  field  has  a 
spiral  angle  of  45°,  that  is,  somewhat  less  radial  than  the 
30°  angle  used  for  Figure  4,  the  draped  fields  should  align 
even  more  with  the  ecliptic  plane  than  illustrated.  On  the 
other  hand,  as  noted  by  McComas  et  al.  [1989],  the  ratio 
of  radial-to-ecliptic  components  does  increase  somewhat 
with  distance  from  the  magnetopause  or  CME  surface. 
But  the  continued  dominance  of  the  ecliptic  component 
out  to  the  shock  leads  us  to  conclude  that  the  importance 
of  out-of-the-ecliptic  draping  mu.st  be  at  most  secondary, 
at  least  near  equinox,  compared  to  the  Russell-McPherron 
effect,  especially  since  only  the  latter  can  contribute  to  the 
semiannual  variation  of  storm  occurrence. 

Discussion 

The  preceding  section  demoastrates  how  strong  south¬ 
ward  IMF  can  arise  from  compressed  ecliptic  fields 
in  the  postshock  flow  of  CMEs  during  equinox.  By 
themselves,  these  fields  could  be  responsible  for  the 
pronounced  semiannual  variation  in  great  geomaignetic 
storms.  However,  they  may  also  serve  as  stepping  stones 
to  the  generation  of  great  storms  during  equinox  by 
priming  Earth’s  ring  current  so  that  additional  southward 
fields  in  the  driver  gas  carry  D^i  over  the  250  nT 
threshold.  Priming  by  postshock  fields  occurs  in  four 
of  the  ten  less  intense  storms  {D,i  <  100  nT)  analyzed 
by  Tsurutani  et  al.  [1988]  with  solar  wind  data,  and  the 
compound  nature  of  the  £).,(  profile  for  two-thirds  of 
the  great  storms  listed  by  Clivcr  et  al.  [1992]  suggests 
that  priming  may  be  even  more  important  for  the  largest 
storms.  Thus,  the  postshnek  Russell-McPherron  effect 
may  serve  as  the  added  factor  that  brings  storms  up  to 
great  .storm  status. 

Yet  whether  or  not  the  postshock  Ru.ssell-McPherron 
effect  can  account  for  the  whole  of  the  pronounced 
semiannual  variation  in  great  storm  occurrence  is  not 
dear.  Tsurutani  et  al.  [1992]  find  that  two  of  their  five 
great  storms  are  caused  by  driver  gas  alone.  If  driver 
gases  alone  are  responsible  for  a  large  fraction  of  all  great 
storms,  as  suggested  by  the  limited  study,  then  it  seems 
that  the  occurrence  frequency  of  driver  gas  storms  must 
necessarily  also  have  a  semiannual  variation,  in  order 
to  account  for  the  fact  that  essentially  no  great  sto.ms 
occur  at  the  solstices.  The  only  other  possibility  is  that 
the  magnetosphere  itself  loses  its  coupling  efficiency  at 
solstices  [e.g.,  Crooker  and  Siscoc,  1986b]. 

A  way  in  which  driver  gas  storms  may  contribute 
to  the  semiannual  variation  is  by  retention  of  some  of 
the  ordered  coronal  fields  in  CMEs.  Hoeksema  and 
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Zhao  [1992]  demonstrate  that  field  orientation  retention 
appears  to  hold  for  four  out  of  five  cases  of  driver-gas- 
driven  storms  analyzed  by  Tsurutani  et  al.  [1988,  1992). 
Theoretically,  if  CMEs  emerge  from  helmet  streamers  at 
sector  boundaries  with  high  inclination,  as  pictured  in 
Figure  3,  and  the  helmet  fields  form  the  leading  edge 
of  the  driver  gas  [e.g.,  Hundhausen,  1988],  then  the 
leading  edge  will  have  roughly  the  same  orientation  as 
the  overlying  draped  fields  and,  like  them,  be  subject 
to  the  Russell-McPherron  effect.  The  fact  that  sector 
boundciries  tend  to  have  high  inclination  near  solar 
maximum  [Hoeksema  and  Scherrer,  1986],  when  CME 
occurrence  is  highest  [e.g.,  Webb,  1991],  gives  weight 
to  the  argument.  Unfortunately,  the  limited  amount  of 
interplanetary  data  available  for  great  storms  hampers 
direct  testing  of  this  hypothesis. 

Conclusions 

1.  The  semiannual  variation  of  great  geomagnetic 
storm  occurrence  and  of  monthly  avarages  of  Dgt  in  years 
of  great  storms  Is  remarkably  pronounced.  This  implies 
that  the  semiannual  variation  found  by  averaging  indices 
over  many  years  is  rujt  a  statistical  effect  superposed  upon 
random  storm  occurrence  but  primarily  the  effect  of  the 
storms  themselves. 

2.  Compression  and  draping  within  the  ecliptic  plane 
in  postshock  flow  preceding  CMEs  greatly  strengthens 
the  Russell-McPherron  effect  there,  which  can  account 
at  least  in  part  for  the  semiannual  variation  of  great 
magnetic  storm  occurrence. 

3.  Whether  or  not  geomagnetic  storms  caused  by 
southward  magnetic  fields  within  the  CME  driver  gas 
contribute  to  the  semiannual  variation  remains  to  be 
tested. 
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